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Abstract. The aim of this work is to derive sharp quantitative estimates of 
the qualitative convergence results developed in 28 for regularized full- and 


partial-cloaks via the transformation-optics approach. Let r 0 be a compact set 
in R 3 and Tj be a ^-neighborhood of To for S £ R+. U represents the virtual 
domain used for the blow-up construction. By incorporating suitably designed 
lossy layers, it is shown that if the generating set To is a generic curve, then one 
would have an approximate full-cloak within S 2 to the perfect full-cloak; whereas 
if r 0 is the closure of an open subset on a flat surface, then one would have an 
approximate partial-cloak within S to its perfect counterpart. The estimates 
derived are independent of the contents being cloaked; that is, the cloaking 
devices are capable of nearly cloaking an arbitrary content. Furthermore, as 
a significant byproduct, our argument allows the relaxation of the convexity 
requirement on r 0 in 28 , which is critical for the Mosco convergence argument 
therein. 
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1. Introduction 


This work concerns the invisibility cloaking for time-harmonic scalar waves gov¬ 
erned by the Helmholtz system via the approach of transformation optics. A 
proposal for cloaking for electrostatics using the invariance properties of the con¬ 
ductivity equation was pioneered in 21 22 . Blueprints for making objects in¬ 


visible to electromagnetic (EM) waves were proposed in two articles in Science 
in 2006 27\ 38 . The article by Pendry et al uses the same transformation used 
in 21,22 while the work of Leonhardt uses a conformal mapping in two dimen¬ 


sions. The method based on the invariance properties of the equations modelling 
the wave phenomenon has been named transformation optics. There have been 
several other proposals for cloaking 12,361. The transformation optics has received 
a lot of attention in the scientific community and the popular press because of the 
generality of the method; see (TTJ[T9, 20, 33],|42] for comprehensive surveys on the 
theoretical and experimental progress. 
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Using the transformation optics, the construction of an ideal cloak is based on 
blowing up a point in the virtual space which splits a ‘hole’ in the physical space to 
form the cloaked region, whereas the ambient space is ‘compressed’ via the push- 
forward to form the cloaking region. In a similar fashion cloaking devices based on 
blowing up a crack (namely, a curve in M 3 ) or a screen (namely, a flat surface in 
M 3 ) were, respectively, considered in 13,18 and 29 , resulting in the so-called EM 
wormhole and carpet-cloak respectively. The blow-up-a,-point (respectively, -crack 
or -screen) construction yields singular cloaking materials, namely, the material 
parameters violate the physical regular conditions. The singular media present 
a great challenge for both theoretical analysis and practical fabrications. In or¬ 
der to avoid the singular structures, several regularized constructions have been 
developed in the literature. In 15 17,40 , a truncation of singularities has been 
introduced. In 25, 26,31 , the blow-up-a-point transformation in 22,27,38 was 


regularized to become the ‘blow-up-a-small-region’ transformation. Nevertheless, 
it is pointed out in 24 that the truncation-of-singularities construction and the 
blow-up-a-small-region construction are equivalent to each other. Instead of the 
ideal invisibility, one would consider approximate/near invisibility for a regular¬ 
ized construction; that is, one intends to make the corresponding wave scattering 
effect due to a regularized cloaking device as small as possible depending on an 
asymptotically small regularization parameter 5 G R+. 

Due to its practical importance, the regularized cloaks have been extensively 
investigated in the literature; see 130 [25) [26) [32) 


34 and the references 


therein. In a recent work 28 , a mathematical framework on constructing reg¬ 
ularized full- and partial-cloaks for time-harmonic scalar waves governed by the 
Helmholtz system was developed, extending the existing study in the literature to 
the most general case. There are three major ingredients of the general cloaking 
scheme in 28 . The cloaking medium is obtained by blowing up a virtual domain 
T s , which is a ^-neighborhood of a compact convex generating set T 0 . Second, a 
compatible lossy layer should be employed right between the cloaked and cloak¬ 
ing regions, which is indispensable of enabling the corresponding construction to 
nearly cloak an arbitrary optical contents. Third, based on some qualitative con¬ 
vergence analysis, it is shown that if the generating set To has a zero capacity, then 
one would have an approximate full-cloak; whereas if the generating set To has a 
non-zero capacity, then one would have an approximate partial-cloak. Here, by a 
full-cloak, we mean that the invisibility holds for incident waves coming from every 
possible direction and observation made in every possible angle, and otherwise the 
cloaking device is called a partial-cloak. 

Clearly, quantitatively estimating the convergence results in 


28 would yield 


a much more accurate characterization of the regularized cloaking constructions. 
Most of the available results in the literature mentioned earlier on quantitatively 
characterizing the degree of approximation to the ideal cloak are mainly concerned 
with the case that the generating set To is a single point. For the other cases, the 
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corresponding mathematical study in deriving sharp quantitative estimates were 
fraught with significant technical difficulties and were lacked in the literature. This 
work aims at dealing with this challenging issue. Specifically, by letting T 0 and Tj, 
5 G M+, respectively denote the generating set and virtual domain, we show that 
if T 0 is a generic curve, then one would have an approximate full-cloak within 5 2 
to the perfect full-cloak; whereas if To is the closure of an open subset on a flat 
surface, then one would have an approximate partial-cloak within 5 to its perfect 
counterpart. The estimates derived are independent of the contents being cloaked; 
that is, the cloaking devices are capable of nearly cloaking arbitrary contents. 
We make essential use of variational and layer potential techniques for our study. 
Some of the results derived in the process are of significant theoretical and practical 
interests for their own sake. Furthermore, as a significant byproduct, our argument 
allows the relaxation of the convexity requirement on T 0 in 
for the Mosco convergence argument therein. 

The rest of the paper is organized as follows. In the next section, preliminaries 
on acoustic scattering and transformation-optics cloaking are presented. In section 
3, we introduce the layer potentials and briefly discuss their mapping properties. 
Section 4 is devoted to showing the sharp quantitative estimates of the regularized 
full-cloak when the generating set is a generic curve. In section 5, sharp estimate 
on the regularized partial-cloak when the generating set is a flat subset. The main 


which is critical 


results in this paper are given in Theorems 4.1 and 5.1 


2. Wave scattering and transformation-optics cloaking 

Let q G L°°(M 3 ) be a (complex valued) scalar function and a = (cP 7 ) G L°°(M 3 ) 3x3 
be a real symmetric-positive-definite matrix-valued function. It is assumed that 
there exists 0 < A < 1 such that 


Jh/(x) > A, 9g(x) > 0; AI 3x3 < cr(x) < A 1 1 3x3 for a.e. x G M, (2.1) 

where I 3x3 is the 3x3 identity matrix. Physically, q is the modulus and a -1 is the 
density tensor of an acoustic medium in M 3 . It is always assumed that the acoustic 
inhomogeneity is compactly supported, i.e., there is a bounded domain C M 3 
such that q = 1 and cr = I 3x3 in M 3 \fi. We write {fi; a, q] to signify an acoustic 


medium as described above. (2.1) is the regular condition for the acoustic medium. 
In scattering theory, time-harmonic plane waves of the form u*(x;o;,d) := e ia;x ' d 
impinge on the acoustic medium {fi; a, q} from the direction d G S 2 := {x G 
M 3 ; |x| = 1}. Here, c <j G M+ denotes the wave number of the time-harmonic wave. 
The corresponding total wave field w(x;o;,d) G C, x := (xj) 3 =1 G M 3 , is governed 


by the following Helmholtz system 

3 ° 7 d 



( cr J (x)—u(x; d,(n)J + ca g(x)w(x; d, ca) =0, x G M , 

:= u(x) — wTx) satisfies the Sonnnerfeld radiation condition. 


( 2 . 2 ) 
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IUU 


= 0 , 


(2.3) 


By the Sommerfeld radiation condition, we mean that the scattered wave n + (x) 
satisifies 

lim Ix|f a " +(X) 

|x|A+oo 1 1 V 0|x| 

which holds uniformly in x := x/|x| G S> 2 . The radiation condition characterizes 
the decaying property of the scattered wave at infinity. The forward scattering 
problem ( 2.2[ ) has a unique solution u G Hl oc (M. 3 ) (cf. 135] ) and moreover, as 
|x| —* +oo (cf. [121) 


u(x;d,cu) = e iwx ' d + 


Ju\x\ 


~U r 


(x, d) + O 


(2.4) 


which holds uniformly in x G § 2 . Moo(x, d) is known as the scattering amplitude , 
and it is real-analytic with resect to x and d, which are respectively referred to 
as the incident direction and observation angle. Clearly, u 00 depends on the inci¬ 
dent wave field e luJX ' d as well as the acoustic medium {12; a, q}. In what follows, 
we shall also write d; {12; a, q}) to signify such dependences, noting that we 

have dropped the dependence on u>. d; {12; a, q}) for all (x, d) G § 2 x § 2 

carries the optical information of the acoustic inhomogeneity {12; a, q}. An im¬ 
portant inverse scattering problem arising in practical applications is to recover 
the medium {12; a, q} by knowledge of , u 0O (x, d). This inverse problem is of fun¬ 
damental importance to many areas of science and technology, such as radar and 
sonar, geophysical exploration, non-destructive testing, and medical imaging to 


name just a few; see 12,23,42 and the references therein. In this context, an ideal 


or perfect invisibility cloak is introduced as follows. 


Definition 2.1. Let 12 e and I2j be bounded domains such that 12* d 12 e . f2 e \f2j 
and f2j signify, respectively, the cloaking region and the cloaked region. Let A 0 and 
Aj be two subsets of § 2 . {f2 e \I2j; a c , q c } is said to be an (ideal/perfect) invisibility 
cloak for the region 12* if 

Uoc (x, d;{I2 e ; cr e ,g e }) = 0 for x G A„, d G A;, (2.5) 


where the extended medium is given as 


{I2 e , cr e , q e { 


Oa, Qa ill 

qc in f2 e \f2 i 7 


( 2 . 6 ) 


with {12*; a a , q a } denoting a target medium located inside f2j. If A 0 = A * = S 2 , then 
it is called a full cloak , otherwise it is called a partial cloak with limited apertures 
A 0 of observation angles, and A* of impinging angles. 


By Definition 2.1, one has that the cloaking layer {f2 e \fb; a c , q c } makes the tar¬ 


get medium {fb; a a , q a } invisible to the exterior scattering measurements when the 
detecting waves come from the aperture Aj and the observations are made in the 
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aperture A c . From a practical viewpoint, the cloaking device should be indepen¬ 
dent of the target medium <j a , q a }. Indeed, throughout the present study, the 

target medium {hip a a , q a } could be arbitrary (but regular). As discussed earlier in 
the introduction, the cloaking medium {fi e \hh; a c , q c } for an ideal cloak are usually 
singular; namely, the regular condition (2.1) is violated. Next, we introduce the 
definition of an approximate cloak by making use of regularized cloaking medium 
{n e \rh; erf, gf }, which is the main theme of the current work. Here, 5 G M + signi¬ 
fies a regularization parameter, and {h2 e \hh; erf, gf} approaches {0 e \0*; cr c , g c } as 
5 —» +0 in a certain sense. 


Definition 2.2. Let H e ,h2j and A 0 ,Aj be as introduced in Definitio n |2.1| Let 
{D e ;crf,gf} be as defined in (2.6) with {D e \Dj; a c , q c } replaced by {D e \Qj; erf, gf}. 
Here, {D e \fh; erf, gf } is as described earlier, and it is called an approximate cloak 
with incident aperture A, ; and observation aperture A 0 for the region Q, if 

Woo (x, d; {D e ;crf,gf}) = £(S) for x G A 0 , d G A*, (2.7) 

where £(S) is a positive function satisfying £(S) —> 0 as 5 —* +0. 


Clearly, in Definition 2.2, £(S) characterizes the accuracy of the regularized ap¬ 


proximate cloak. In our subsequent study, we sometime do not explicitly specify 
whether a cloak is ideal or approximate, and it should be clear from the con¬ 
text. In order to construct an invisibility cloak through the transformation-optics 
approach, a critical ingredient is the following transformation properties of the 
optical parameters (cf. [16,28]). 


Lemma 2.1. Let L2 and D be two Lipschitz domains, and suppose there exists a 
bi-Lipschitz and orientation-preserving mapping 

x = F(x) : Q —> Q. 

Let (D; a, q) be an acoustic medium supported in D. Define the push-forwarded 
medium as follows, 

{D; a, q} := F*{D; a, q}, (2.8) 

where 


cr(x) = F*cr(x) := (fif(x) • cr(x) ■ F>F(x) T ) |det(F>F(x)) | 
g(x) = F*g(x) := g(x)|det(-DF(x ))| -1 


-i 


c =F-!(x) 


(2.9) 


x=F- 1 (x) 


where DF denotes the Jacobian matrix of F. Then u G i/ 1 (D) solves the Helmholtz 
equation 

V • (crVu) + ui 2 qu = 0 in D, 

if and only if the pull-back field u = (F” 1 )*!! u o F G i7 1 (D) solves 

V ■ (cfVu) + l o 2 qu = 0 in LI. 
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We have made use of V and V to distinguish the differentiations respectively in 
x- and Si-coordinates. 

As a direct consequence, one has that if F : is a bi-Lipschitz mapping 

such that F\gtt = Id, then 


Ur 


>(x, d; {fi; a, q}) = Uoo(x, d; F*{H; a, q}). 


( 2 . 10 ) 


As a general description, in our subsequent study, 5 G M + always denotes the 
small regularization parameter. Let O e , ffj be the bounded domains in Defini- 
and 


tions 


2.1 


2.2 


We shall also need D $/ 2 d D$ d hi* and hi* d hi; d hl e , which 
shall be explicitly specihed in our subsequent study. D$/ 2 and D, 5 are the virtual 
domains, and actually plays the role of in our earlier discussion. Throughout 
the rest of the paper, we assume that there exists an orientation-preserving and 
bi-Lipschitz mapping F$ : Q e —> such that 

Fs(Ds/2 ) = &i, Fs(Ds\D$/ 2 ) = F$(D e \ D$ ) = Q e \Q t and Fs\gci = Id. 

( 2 . 11 ) 

We refer to 128 for the discussion on the so-called assembled-by-components strat¬ 
egy of explicitly constructing such a transformation for generic D$ and Vt e . In our 
study, an approximate cloaking structure {h 2 e ; cq q} hi the physical space is ob¬ 
tained by pushing forward a virtual structure {fl e ; <7s,qs}, 

/, 1 in Q e \ D, 5 , 


{fl e -,cr s ,qs} = (F s )*(O e ;c 7 ,g) = { a u qi in D 5 \D 5/2 , (2.12) 

a a , q a in D s/ 2 , 

where a a , q a are arbitrary but regular. Set 

{fl e ; (T, q} = (F s )*{n e -, a s , qs} and {D 5 / 2 -, o a , q a } = (F^ 1 )^^; a' a , q a }, 

which denote, respectively, the physical cloaking structure and the physical object 
being cloaked. Clearly, {hh; a' a , q' a } is also arbitrary but regular. The design of eq 
and qi in the lossy layer Ds \ Ds/ 2 is critical and will be explicitly given in what 


follows at the appropriate place. By the transformation acoustics in Lemma 2.1 


it is straightforward to show that in order to quantify a cloaking device {fl e ; a, q } 
in the physical space, it suffices to derive the approximate cloaking effect to the 
virtual structure {Q e ] a s , q s }. In other words, it is sufficient to show that the 
scattering amplitude u ^ (x, d) to the following system in the virtual space 


Aus + tv 2 u$ = 0 in M 3 \ Ds, 

V • a 6 Vus + u 2 q 5 = 0 in Ds, 

u + := us — u l satisfies the radiation condition (2.3), 


(2.13) 


verifies (2.7). Hence, throughout the rest of our study, we shall focus on quantifying 


the scattering of the virtual system (2.13). 
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3. Boundary layer potentials 


In this section, we introduce the boundary layer potential operators for our 
subsequent use, and we refer to |8[|l2,35^37,4l] for more relevant discussions. Let 

l Kn r nn ninnomnurn cn mrirvn i /~\ /\ 1 /.i" mi rr\ Kt 


be the fundamental solution to A + or given by 


GUx) = 


0 iuj |x| 


47t|x| 


(3.1) 


For lu = 0 the corresponding fundamental solution is denoted by G'(x). Let D be 
a Lipshitz domain. Denote by Sp the single layer potential operator 

«S£[0](x) : = f ^u( x — y)0(y ) dcr y x6l 3 \ad. (3.2) 

JdD 

It is well-known that the single layer potential S satisfies the trace formula 
|;SSM| ± = (±b + (K^y) [0] on dD. 

Here and throughout the rest of the paper, we make use of the following notation. 
For a function u defined on M. 3 \dD , we denote 


(3.3) 


and 


du 

dv 


«L(x) = lim w(x ± r ■ p(x)), x G dD, 

t—H-0 


(x) = lim (V x u(x ± r • i/(x)), z/(x)), x G dD, 

T— s>+0 


if the limits exist, where i/(x) = v x is the unit outward normal vector to dD at 
x G dD. In (3.3), (/Q))* is the L 2 -adjoint of /Cp and 


: = P- v - 


9G U (x — y) 


0(y) da y , x G dD 


JdD du(y) 

where p.v. signifies the Cauchy principle value. We also define the double layer 
potential operator by 

r dG u (x — y) 

IdD 


^d[0 ]( x ) : = 


dv ,, 


-0(y) do y xgB 3 \ dD 


(3.4) 


and we have the following trace formula 




± 


= Tx/ + /C‘ 


B 


on dD. 


(3-5) 


By using the layer potential operators introduced above, we make use of the 

(3.6) 


following integral representation to the solution of (2.13) outside of D§ 

l 3 \D~s, 


us = u l + S' 


D s [ 


111 
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where it is easily verified that the density function <f> satisfies 


i + ^ Ds r)w = 8{us ~ u ' 


du 


dD s 


(3.7) 


4. Regularized full-cloak 

In this section, we consider the regularized cloaking construction by taking the 
generating set to be a curve. Let T o denote a smooth simple and non-closed curve 
in M 3 with two endpoints, denoted by P 0 and Q 0 , respectively. Let r £ M + . Next, 
by using r 0 as a generating set, we shall construct a simply connected set D r in 
M 3 for our cloaking study; see Fig. [I] for a schematic illustration. Denote by 7V(x) 
the normal plane of the curve r 0 at x £ r 0 . We note that N(Pq) and N(Qq) are, 
respectively, defined by the left and right limits along r 0 . For any x £ r 0 , we let 
J^ r (x) denote the disk lying on 7V(x), centered at x and of radius r. It is assumed 
that there exists ro £ R+ such that when r < ro, ^(x) intersects To only at x. 
We start with a thin structure D{. given by 

D f r :=y r {x) xr„(x), x £ r 0 , (4.1) 

where we identify r 0 with its parametric representation r 0 (x). Clearly, the facade 
of D/, denoted by S'/ and parallel to r 0 , is given by 

S'/ := (x + r • n(x); x £ r 0 , n(x) £ N(x) D § 2 }, (4.2) 

and the two end-surfaces of Df are the two disks S^ t {Pq) and 5? r (Q o)- Let 
and D b be two simply connected sets with dD^ o = S*/ o U ^ ro (P 0 ) and dD b o = 
S b Q U 5? ro {Qo)- It is assumed that S ro := S'/ U S b 0 U S'® is a C' 3 -smooth boundary 
of the domain D ro := D/ o U D£ o U D b rQ . For 0 < r < r 0 , we set 

D a r ■= -(D a ro -P 0 ) + P 0 = (- • (x-P 0 ) + ^o;x £ D a \ , 

r 0 0 [r 0 ) 

and similarly, D b := r/r 0 ■ (D b o —Qo) + Qo- Let S'/ and S b , respectively, denote the 
boundaries of £)/ and D b excluding and 5^ b . Now, we set D r := D“UD/uHj, 
and S r := S'/ U S b U S'/ = dD r . 

According to our earlier description, it is obvious that for 0 < r < ro, D r is a 
simply connected set with C 3 -smooth boundary S r , and D ri d D r2 if 0 < r\ < 
t 2 < r 0 . Moreover, D r degenerates to T 0 if one takes r = 0. In order to ease 
the exposition, we assume that r 0 = 1. In what follows, we let 5 £ M + be the 
asymptotically small regularization parameter and let 

S s :=S f s US b s USt, (4.3) 

denote the boundary surface of the virtual domain Ds used for the blowup con¬ 
struction. In order to ease the exposition, we drop the dependence on r if one 
takes r = 1. For example, D and S denote, respectively, D r and S r with r = 1. It 
is emphasized that in all our subsequent argument, D can always be replaced by 
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FIGURE 1. Schematic illustration of the domain D r for the regular¬ 
ized full-cloak. 


D tq with 0 < r 0 < r 0 being a fixed number. Hence, we indeed shall not lose any 
generality of our study by assuming that r 0 = 1. Finally, we would like to note 
that a particular case is to take r 0 to be a straight line-segment and, S% and S “ 
to be two semi-spheres of radius r and centered at P 0 and Q 0 respectively, which 
is exactly the one considered in 28 . Hence, the geometry in our study is much 


more general than that considered in 28 ; see also Remark 4.4 in what follows for 
more relevant discussions. 

We are now in a position to present the main theorem on the approximate cloak 
constructed by using D§ described above as the virtual domain. 


Theorem 4.1. Let D$ be as described above with c)D d = S $ defined in (4-ty - Let 
us be the solution to (2.13) with {12 e ; as, qs} defined in (2.12), and {Ds\Ds/ 2 \ &i, Qi} 
given in (4.11). Let u^fk, d) denote the scattering amplitude of us- Then there 
exists 5 0 G M + such that when S < 5 0 , 

|»i(i,d)|<CM« ! , (4.4) 

where C(oj) is a positive constant depending on u> and D, but independent of a a , 
q a and x, d. 

Remark 4.1. Following our discussion at the end of Section 2, one can immediately 
infer by Theorem |4.1 that the push-forwarded structure 

{fl e ; a , q} = (F 5 )*{Q] a s , qs} 


with Fs defined by ( 2 . 11 ) and { 12 e ; as, qs} defined in the theorem, produces an 
approximate cloaking device within h 2 -accuracy to the ideal cloak. Indeed, if one 
lets Moo denote the scattering amplitude corresponding to the cloaking structure, 
then there holds 

|-Uoo(x, d)| < C(u)S 2 (4.5) 

where C(oj) is independent of o' a , q' a and x, d. 
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Remark 4.2. It is remarked that our estimate derived in Theorem |4.1| is sharp and 
this can be vended by the numerical examples presented in 


28 


The subsequent three subsections are devoted to the proof of Theorem 4T For 
our later use, we drst derive some critical lemmas. In what follows, we let z denote 
the space variable on T 0 and for every y G Df we dedne a new variable z y G T 0 
which is the projection of y onto To- Meanwhile, if y belongs to D “ (respectively 
D b r ) then z y is dedned to be Po (respectively Q o). Let t denote the arc-length 
parameter of T 0 and 6 , which ranges from 0 to 2i r, be the angle of the point on 
N(t ) with respect to the centeral point x(t) G T 0 . Moreover, we assume that if 
6 = 0 , then the corresponding points are those lying on the line that connects 
T 0 (t) to T i(t), where Tx is dedned to be 

Ti := {x + ni(x); x G T 0 , ni(x) G IV(x) D § 2 — a dxed vector for a given x}. 


Next we introduce a blowup transformation which maps y G D§ to y G D as 


follows 






4f(y) = y: 

= ^(y-Zy) +Z y , 

y e D f s , 

(4.6) 

whereas 

4l(y) = y 

._ J r A + p °. 

' \ *y + Q„, 

y e Df, 
y e D b s . 

(4.7) 


Then, we can show that 


Lemma 4.1. Let A be the transformation introduced in (4.6) and (4.7) which 
maps the region D$ to D. Let B(y) be the corresponding Jacobi matrix of A( y) 
given by 

B( y) = V y H(y), y eWs. (4.8) 

Then we have Biy) 1 = B(y) and 


B{y)v y = -Uy, y G dD Sl 

where v y stands for the unit outward normal vector to dD$. 
Proof. It is easily seen that if y G D$ U D$, then 


(4.9) 


B (y) = ^I 3 X3, y G Dg U D b . 

Next, we let y G and note that y — z y = S(y — z y ), namely 

^(y) = y = Jy- (\~ l ) z y 

Hence it follows 

B( y) = V y H(y) = ^I 3x3 - (J - l)z' y (t)z' y (t) T , (4.10) 
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where the superscript T denotes the transpose of a vector or matrix. Then it is 
straightforward to see that 


B{y)v y = ^Uy, y G dD s . 


The proof is complete. 


□ 


Next, we present the crucial design of the lossy layer. Define the material 
parameters oy and q$ in the lossy layer D$ \ D$/ 2 to be 

cr 5 (x) = cr z (x) := 7 -5~ 2 (x), q s (x) = ^(x) := a+i/3 for x G D 5 \D 5/2 , (4.11) 


where 7 , a and j3 are all positive real numbers. It is easily seen from (4.10) that 
B is a symmetric positive definite matrix and hence 07“ 1 is a well-defined regular 
material tensor. 


4.1. Asymptotic expansions. In order to tackle the integral equation (3.7), 
we shall first derive some crucial asymptotic expansions. Henceforth, we denote 
0(y) := 0(y) for y = A(y), y e D$ and y G dD. The same notation shall be 
adopted for x and x. First, we note that 


day = 1 sda «’ 


6 S{, 


VdOy, yeSJUS?. 


(4.12) 


For a function u, its gradient in D J can be decomposed into two parts: the tan¬ 
gential derivative with respect to r 0 and the corresponding normal derivative as 
follows 

V„tt(y) = d t u(y)T y + V/w(y), 

where T y is the tangential direction along To. Here, V; signifies the normal de¬ 
rivative. Then for x G M 3 \ D$ with sufficient large |x|, we can expand [</>] in 
z G r 0 as follows 


S 


D S 


X = 


G u (x-y)(p(y)da y 


1 dD s 


1 dDs 


(g u (x - z y) - VG w (x - z y) ■ (y - Z y )^0(y) do y 


+ 


(y 


\T V 72 


'dDs 


>dDs 

G w (x - z y )0(y) da y - 


V 2 Gh,(x - C(y))(y - z y )<f>(y)da y 

VG w (x - z y ) ■ (y - z J/ )0(y) da y 


'dDs 


+ 

■=Ri 


(y - z y ) T v 2 G w (x - C(y))(y - z y )<f>(y)d(r y 


dD s 

F R 2 




(4.13) 
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where £(y) = rjy + (1 — r])z y G Ds for some rj G (0,1), and the superscript T 
signifies the matrix transpose. The term i? 3 in (4.13) is a remainder term from 
the Taylor series expansion and it verifies the following estimate 


m = s 2 


i 


IdDg 


(y-z») v Gu(x - C(y))(y ~ Zy)0(y)da p <6 rj|M| H -j/j( 8£)) . 

(4.14) 

We shall also need the expansion of the incident plane wave u l in z G T 0 , and there 
holds 


«*(y) = u \ z y) + Vu \ z y) ' (y ~ z y)+Y ^u t (z y )(y - z y )°, (4.15) 

l«l=2 

where the multi-index a = (or, a 2 ,a 3 ) and <9" = <9" 1 <9,“ 2 . Since for y G dD$, 
u y = Uy, one further has 

OO OO 

v v ■ v X(y) = v v • v Y d y ul ( z v )(y - z y) a = Y 5H ~ ld y ul ( z y) u v ■ v (y - z v) a - 

|a|=0 |a| = l 

The following lemma is of critical importance for our subsequent analysis. 


Lemma 4.2. Let (j) be the solution to (3.1) and 0(x) = 0(x) for x G dD$ and 
x G dD. There hold the following results 

mifxl = I + °^ 2 UWH-z/HdD)), x G S f s , 

\ (K:s C ym±) + o(s\m H -s /HdD) ) xe ^, C G{ a ,&>, 

(4.16) 

and 

oo, UUxl = f X G S f s , 

l S(S S c + 5^ / )[0](x) + 0(5 2 \\(j)\\ H -3/2( aD )), x G SJ, c G {a, b}, 

(4.17) 

where 


:= ^ 


To 


( Z X Z 'T V xf _ ( Z x p ioj\z s -Zy\ 


| z x Z j/P 


\%x | 


r»27r 


and 


and 


0(y) dd y ds t , 
(4.18) 


(/C 5 c) h 


x := 


47T 


[ (k-y,v x )~ 

L |x-y| 3 * (y)d<7 *’ 


s• 


D iui |z 2 —z s J /-27r 


'5/ 


x := -- 


r 0 l z s z y\ Jo 


47r 

The variables z x ,z y are inT 0 


(j)(y)d6 y ds t , <S S c[0](x) := / G(x-y)0(y) da 
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Proof. For x, y G D$, one has 

|x y| — |(x z x ) - (y - z y ) + z x - z y \ = |<f((x - z x ) - (y - z y )) + (z x - z 5 )| 
Hence, we have the following expansion for x G S/, 


|x - y| = |zj - Zy\ + <5((x - z £ ) - (y - Zy), 


Z Z,~ 


| Z X Zy| 


r) +0(<5 2 ). 


Similarly 

(x - y, i/j.) = (zjj - z g , I/*) + <5((x - z*) - (y - z s ), i/j,.), 


|x — y| = \Z 


x z y 


1-1 




ZX Zy 


\Zx Zy | 


, (x - z x ) - (y - Zy)) + 0(S 2 ), 


e^'x-yl = e^ z *- z -l f 1 + io;5((x - z,) - (y - Zy), 


Z,T Zy 


\Z£ Zy | 


))+o(s 2 ). 


With those expansions at hand, we can compute for x G SL 


(Wl(x) = s / s;Mx _ y| 


— ICO 


<x-y,^) ^| x _y| 


l x — y| 


0(y) n-3/2(a_D)) 


=<5 i/ d z s _z yl 1 


(z x Zy, V x ) 


\ z x z y | 


,iu;|zg—zg| 


r* 27 T 


M / 0(y) d9y dst + C?(5 11 011i?— 3/2 (azp)), 


which proves the c ase of (4.16) for x G Sg. In a similar manner, one can prove the 
first case of (4.17) with x G Sg. Next, we note that if x G Sg, c G {a, b}, then 

z* - Zy = 0, for y G S c 5 

and 


|x-y| = <J|x-y|, 


(x - y, u x ) (x - y, u x ) 


l x — y I 


l x — y I 


, x, y G Sg. 


Then by using a similar argument to the proof of the first case of (4.16), we can 
prove (4.16) for x G Sg, c G {a, b}. Finally, the second case in (4.17) can also be 
proved following a similar argument by noting that 

l x — y| = S|x — y|? x ,y g Sg. 

The proof is complete. □ 


4.2. An important result for acoustic scattering. In this subsection, using 
the results obtained in the previous subsection, we present a theorem concerning 
the scattering from a thin scatterer. The asymptotic expansion formula derived 
in the next theorem would find important applications in inverse scattering the¬ 
ory. Indeed, it might be used to devise some novel inverse scattering scheme for 
recovering a thin sound-hard obstacle or a thin acoustic medium with a uniform 
content. Hence, the result is of significant practical and theoretical interests for 
its own sake. To our best knowledge, there is no available result in the literature 
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of this form. Moreover, (4.19) derived in the theorem shall be needed in our sub¬ 
sequent proof of Theorem 4T concerning the approximate cloak with an arbitrary 
content being cloaked. 


dl/y 


Theorem 4.2. Letus be the solution to (2.13). Define 4>(y) := 4>(y) = dus(y ^ 
then there holds for x e M. 3 \D, 


9D S 


/ G u (x- Zy)4>(y)day = 2 / G w (x-Zj|)$(y)d(T S +0(5(11^11^-3/2(5/)+ 1)). 
Jsf Jsf 

(4.19) 

If one assumes that <E>(y) = 0, y 6 dD$, namely D$ is a sound-hard obstacle, then 
there holds for x6 K 3 \h 


(u 5 - u l )(x) =2 5 2 / V;G^(x - z y ) • (y - z y )vy ■ Vu l (z y ) day 
Jsf 

- s 2 / G w (x - z y )Aiu l (z y ) day 

Jsf 

- S 2 f Guf* ~ Z y)Vy ■ Vu\Zy)day + O ) , 

Js a us b 


where for a function u, A t u = Vi ■ Wiu. 


Proof. By (3.7), we see that 


0 (x)= (- + (££,)' 


— 1 r + 

v y - V{u 5 -u l )(y) (x). 

oD s J 


Using the results in Lemma 4.2, along with the fact that (cf. |37,41 


- + (JC S c)* : H~ 3/2 {S C ) -A H~ 3/2 (S c ) is invertible, 


and the expansion of u l in (4.15), we have for x e S'/ that 

0(x) = 2$(x) - 2z/ x ■ V,u'(z*) + 0 ( 5 ( 11 ^ 11 ^- 3 / 2 ( 5 /) + !))■ 
By using Green’s formula 


f*2lT p2 7T 

/ v ■ Viu l (z) dd — / l/dd ■Viu l ( z) = 

'o Jo 


3,1) 


(4.20) 


(4.21) 


Vi-V 1 u i (z)ds = 0, (4.22) 


where B(0,1) is the unit disk. Combining (4.21) and (4.22 


, one can show (4.19) 


Next if <E»(y) = 4>(y) = 0, it then can be seen from (4.21) that 

Ill’ll H~ 3 / 2 (dD) < G, 
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where C is a positive constant depending only on D. Then by using our earlier 
results in (4.13) and (4.14), one can first show that 


nuj 

* As 


x = 


' dD s 


G u (x - Zy)0(y) day - / VG u (x - z y ) • (y - z y )0(y) da y + 0(5 3 ) 


'dD s 


=5 / G„ 


'sf 


Zy)(p(y)d<7y + 5 2 


G u 


' S a US b 


Zy)0(y) day 


I sf 


V|G aJ (x - Zy) • (y - Zy)(j){y) day ) + 0(5 3 ). 


(4.23) 


By using Lemma 


4.2 


together with (4.15) and (4.23), we have for x e Sg that 

(4.24) 


0(x) = - 2z/ x • Viu l ( z;g) - 2hz/ x ■ V; ^ (9“u*(z i )(x - z 5 ) a 

| a |=2 

+ A6(^ Sf y[Uy-\/ l U i (Zy)](i) + 0(5 2 ). 

Taking the integral of ( |4.24| ) with respect to 6 from 0 to 27 t and using (4.18) and 
(4.22) one can easily obtain 

/*27T 

/ 0(x) dO = — 2tt5Aiu\z s .) + 0(5 2 ), x G S*. 

Jo 

Therefore, 

[ G„,(x - Zy)4>(y) day = -8 [ G w (x - z y )A l u l (z y ) day + 0(5 2 ). (4.25) 


'sf 


' sf 


Finally, by plugging (4.24) and (4.25) into ( |4.23[ ), along with straightforward 
calculations, one can obtain (4.20) and thus complete the proof. □ 


Remark 4.3. Theorem |4.2| gives the asymptotic expansion of the scattered wave 
held from a thin sound-hard obstacle Ds in terms of the asymptotic parameter 
5 G M+. ffere, by a sound-hard obstacle, we mean a scatterer that the wave cannot 
penetrate inside and the wave velocity vanishes on the boundary of the scatterer 
(cf. (l21). The result can also be extended to deriving the asymptotic expansion of 
the scattered wave held if {Ds] as, qs} is a uniform inhomogeneity; that is, both as 
and qs are two hxed constants. Those results would find important applications in 
inverse problems of reconstructing the scatterers; see |3,4] and references therein 
for related studies. Indeed, as is known that the corresponding inverse scattering 
problems are nonlinear and the asymptotic expansion formulas naturally give rise 
to their linearized counterparts. 

4.3. Proof of the main theorem. Let us go back to the proof of Theorem |4.1| by 
continuing with the estimates of R\ and Ro in (4.13) for the thin virtual cloaking 
structure with arbitrary but regular a a and q a in Ds / 2 . In what follows, we let C 
denote a generic positive constant. It may change from one inequality to another 
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inequality in our estimates. Moreover, it may depend on different parameters, but 
it is independent of a a , q a and x and d. We shall also write C(u>) to signify its 
dependence on the frequency uj. 

Define 

*(y) := *(y) = dus(y) 


dv v 


(4.26) 


We first note that by using (4.21) in Lemma 4.2 there holds 

IHlH-3/ 2 (dD) < C(\\$\\ H -3/2( dD ) + 1). 

Next, by taking expansion around z e To and using (4.19) one can show that for 
6 G M+ sufficiently small and |x| sufficiently large 


(4.27) 


\Ri\ = 
<5 

<25 


dD s 


G„(x - z y )<f>(y) da y 


8 2 


I sf 


G w (x - Zy)cj)(y) day + C— \\<j>\\ H -3/2 (dD) 


(4.28) 


Gb,(x — z^)<h(y) day + C5 2 - r(|| ‘S*IIt/- 3 / 2 (sz?) + 1) 


Jsf 

1 " ~ " ~s:2 1 


—^"mII^II^- 372 ^) i—r(IIIIiT- 3 / 2 (sz?) + !)• 


\R*\ = 


The estimate of i ? 2 can also be done using Taylor’s expansions and (4.21) 
VG^X - z y) ■ (y - Z y )0(y) da y 

I ViG w ( x - z y) • (y - Zy)$(y) day 

'sf 


1 dD s 

<2 h 2 


+ / ViG u (-x-Zy) ■ (y - z y )v y -Viu\zy)day + C8 3 (\\$\\ H -3/2, dD) + 1) 

Jsf 


f2 1 ni^n , i\ , 1 


■p-r(||4 ) || j f/-3/2( S /) + 1) + C5 -p-r(||$||i/-3/2( aD ) + 1). 


(4.29) 


Hence, by applying the estimates in ( |4.14[ ), ( |4.28[ ) and ( |4.29[ ) to (j4T3j) and ( J3.6[ ) , 
we readily have 


I Us — U l | < Ct— r||l»|| fl --3/2(5/) + C5 2 — r(|| < h||//-3/2( aD ) + 1) 


(4.30) 


for |x| sufficiently large. Here, it is noted that we have made use of ||^>||jy- 3 / 2 ^/) 
instead of ||$||tf-i/ 2 (s/) in the above estimates and this shall be crucial in our 
subsequent argument. 


We proceed with the estimate of <f> in (4.30). I 11 the sequel, we set 
u(x) := u 5 (H _1 (x)) = wa(x), xeD. 
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Since 


V ■ a s Vu, 5 + q s uJ 2 u 5 = 0 in D s \ D s/2 , 
by change of variables, one directly verifies that for x G D \ D\/ 2 


7 ( j BVx) ■ (B^BV^usiA- 1 ^)) + (a + ip)u 2 u 5 {A-\Z)) 


0 . 


That is, 

7 Au + (a + i/3)co 2 v = 0 in D\Di/ 2 . (4-31) 

In what follows, we shall estimate ||4 ’||h- 3 / 2 (s/) an d ll < ^ll/ 7 - 3 / 2 (s a us i! )> separately. 
We recall that the H~ 3 t 2 (dD )~norm of the function 4>(x) is dehned as follows 


l$ll H-y^ d D) = sup 

H^Hi/ 3 / 2 (8D)< 1 


$(x)<p(x) da~ x 


' dD 


(4.32) 


Moreover, the H 3 / 2 (5' c )-norm of <f> for c G {/, a, 6 } is given as 


|$| 


H~ 3 / 2 (S C ) SU P 


<1 


H 3/2 (S=)- 


$(x)v?(x) d<Ji 


where H'^ 2 (S C ) denotes the set of i/ 3 // 2 (S' c )-fmictions which have zero extensions 
to the whole boundary dD. We refer to (lj 30,39 for more relevant discussions 
the Sobolev spaces. 

We have 


on 


Lemma 4.3. Let <f> be defined in (4-26), where us is the solution to (2.13) with 
the corresponding as and qs given by (2.12) and (4-11). Then there holds 

II^IIh- 3 / 2 (s/) < + V a2 + /3 2 u 2 ^) \\us\\l 2 (d s \d s/2 ), (4.33) 

and 

ll < 4 , llf7- 3 / 2 (S“us i >) < C5 ^ + \J a 2 + (3 2 uj 2 ^J \\us\\ L 2{ Ds \ D5/2 ). (4.34) 


Proof. For any test function -0 G H^ 2 (S C ), c G {/, a, b }, we let <p be the extension 
of if into dD such that 


f ip on S c , 

\ 0 on dD \ S c . 


(4.35) 


We next introduce an auxiliary function w G H 2 (D) such that 


( w = <p and = 0 on dD, 

i \\w\\h 2 {d) < C'||'0||n3/2( 5 c), (4.36) 

[ w = 0 in D 1/2 U (D \ D c ), 
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The existence of w introduced above can be found in Theorem 14.1, together with 
its Addendum in 


39 


. From the construction, it is readily seen that ^ 


= 0 . 


In the sequel, we let (pi be given as 


dD c 


<Pi = 


0 on S c , 

0 on dD c \ S c . 


Then one has 

/ $(x) 0 (x)dcr 5 = / <F(x)v?i(x) dan + £7 


Vx ■ VxU(x) 


= 7 / Ui • (B ^Vxi'(x) ipi(x)dax + £7 


' dD c \S c 

Vx ■ Vxf(x) 


^l(x) do* 


' dD c \S c 


</t(x) clan. 


Next by using Green’s formula, (4.9), (4.36) and (4.31), one can derive 


<5y 

=^7 

=S / y 


< l(x)' 0 (x) da, 

/ "x ■ Vxufx 


I dD c 


^i(x)d<7 5 


. f dw 

u(x) dan 

w(x) dan ~ / y— 

' V 7 



f dv 
!q D c dvn 

/ Amu — Ainu da 

I d c 


(4.37) 


<h(^V / c?T^a; 2 ||u|| L 2 ( DC \^ / 2 )||w|| L 2 mD c /2 ) 

+ 7|K ; IU 2 (l>7D5 : /2 )||A-u;|| L 2( D c\ D c /2) 

<<5(j + y/V + /3 2 CU 2 ) ||u||L 2 (D 7 D- /2 )lk|| i/ 2 (D) 

<Ch( 7 + y/a 2 + /3 2 (n 2 ) ||u||L2 ( D«=\r»j /a) ||^||H3/2 ( sc ) . 


By change of variables in integrals, it is straightforward to verify that 

lL 2 (D/\Df /2 ) = S 1 W U sWmD f s \D f s/2 y 

’ll L 2 {D c \D c 1/2 ) = d 3/2 11 ^<5 11 z. 2 (T>g\z>| /2 ), ce{a,b}. 


|u| 

\v\ 


(4.38) 


^(x)^(x) dan < C\^y + V W u s\\ L 2 {D f\ D f /2) \\^\\ H ^(sf): ( 4 - 39 ) 

and for c G {a, b} 

$(x)0(x)d(Ti < C<T 1/2 ( 7 +y/a 2 + /3 2 o; 2 ) 11it<511 l 2 (d^\d| / 2 ) 11011 n 3 / 2 (gc). (4.40) 


Finally, we have from (4.37) and (4.38) that 


>sf 
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(4.39) and (4.40) clearly imply that 


II^II 3 / 2 (5^) — c(^y 4“ \/cr 2 4" d 2 ^ ^ II m <5|Il 2 (dJ\d^ 2 )’ 

and for c E {a, b} 

ll^lli?- 3 / 2 (s c ) < C5 1/2 (7 4- \/a 2 + fd 2 u 2 ^ ||m5||l 2 (d^\dj /2 ), 


which further imply (4.33) and (4.34). 
The proof is complete. 


□ 


We next derive the estimate of \\us\\l 2 (d s \d s/2 ) that is needed in (4.33) and (4.34) 
in Lemma [4.31 


Lemma 4.4. Let u$ be the solution to (2.13) , where the corresponding as and qs 
are given by (2.12) and (4-H)- Then there holds 

IMILd,\d„j) < ■^|i»o 0 |'+Ci5^5||'I|| i ,-3/j (i 5/ l +CT 2 — \\i\\ H -W{SD)+cr-—, 

(4.41) 

where C is independent of u. 


Proof. By (2.4) and the Sommerfield radiation condition 


d(us — u l ) 


dr 


= iu(us — u l ) + 0( |x| ), as r = |x| —> 00, 


it is straightforward to show that 

9 lim / 


R—t 00 


' 0 B R 


d(u s - u l ) 

(u s ~ u l ) 


dr 


= uj lim 

R—> 00 


' 0 B R 


\u s - u l I 2 , 


(4.42) 


where Br denotes a central ball of radius R. Next, noting that 

A u l + a fu 1 = 0 in R 3 , 


we multiply both sides of the equation (2.13) in R 3 \ D$ with us — u\ and in D$ 
with us, and integrate by parts to have 


lim 

R —>00 


d(u 5 - u l ) _• 

(u s - 'id) - 


1 dB R dr 

d(u s - u*) 


lM. 3 \D s 


1 dD s 


dv 


(us — u*)|+ 4- / u-asVus 
J dDs 


\\ 7 (u s - u l )\ 2 + u 2 / _ | u s - u l \ 2 

Jr 3 \Ds 

Us - ■ a s Vus 


'As 


+ u 2 qs\u 5 \ 2 = 0 . 
Jd s 
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By taking the imaginary parts of both sides of the last equation and using the 

(4.43) 


definitions of as and q$ in (2.12) and (4.11), we arrive at 

kl 2 




R—>oo 


8 ^- u \u s -v‘)=-^ 


'OBr 


dr 


3 


' D S \D S /2 

du 5 


'dD s 


dv 


du 


u'-S —(u s -u l ) | + . 


'dD s 


dv 


We next estimate the last two integral terms in the RHS of (4.43). By the definition 
of $ in (4.26), one has 


r du s 
JdDs dv 

u l 

= 

[ $v* 

< C5 

[ rT(z)$(y )dxTy 

+ 


JdDs 


Jsf 


(4.44) 


<C5\m H - 3/2{sf) + C5 2 \\§\\ H - 3 /2( dD y 
Similarly, by (|3.6|), (|4.17l) and (|4 . 27[) , one has 


[ d -f{us-u l ) | + 
I dDs dv 
c2 



f du 1 r . 

— 

L 


(4.45) 


<C5 2 \\(j)\\H-3/2(dD) < C5 2 (\\$\\ H -3/2( dD ) + 1 ). 

Finally, by combining (|4.42|), (4.43), (|4.44[) and (4.45), we have 


/3uj I _ kl < cj ^lini j \us — u l \ + C5 114*11 H~ 3 / 2 (sf) 4“ C5 (|| < h||//-3/2(g£ ) ) + 


' D s \D 5 /2 


' 9Br 


which together with the fact that 


lim 

R— xx) 


k — u l \ 2 < C\u 0 


9Br 


readily yields (4.41). 

The proof is complete. 


□ 


Proof of Theorem f.l The key ingredient of the proof is to derive the following 
two estimates 


ll^k-3/2(s/) < C(u)5 and \\$\\ H -3/2 {dD) < C(u). 


(4.46) 


From (4.34) one can obtain 

12 


\\m ~ 3 / 2 (dD) — cm\ 2 H -3/2(5/) + kk- 3 / 2 (S“US i >) ) < C{uf)5 \\'98\\l2(Ds\D s/2 )- 

(4.47) 


By (2.4) and (4.30), we have 


Uoo (±,d)\<C 5||$|| 

R-3/2(5/) + <5 (\\®\\h- 3 / 2 {9D) + 1) )> 


and therefore 


^ockd)! 2 < cm\*r H -*, 2 {S f) + s3 m 2 H- 3 / 2 (9D) + t 


(4.48) 

(4.49) 




































REGULARIZED FULL- AND PARTIAL-CLOAKS 


21 


Then by plugging (4.41) and (4.49) into (4.47), one can derive 
12 


\\®\\H-3/2(dD) (JUool + ^II < ^ ) ||r- 3/2( S /) + h ||$|| H -3/2( 9D) + h 

<C(u) (<5||i>||^ -3/2(5/) + ^ 2 ||^||^-3/2 (aD ) + ||$||i/-3/2(5/) + £||$||//-3/2(9 D ) + 

<C'(o/)( ( 5||<h||^ 

-3/2 (dD) + 11*11* -3/2(9.D) + ^)) 

which readily implies by choosing 5 G M + sufficiently small that 

11^11/5-3/2(90) < C{u). (4.50) 


Next by (4.33) and (4.41), along with the use of (4.50), there holds 


I $11/5-3/2(5/) <C(u) Muool 2 + 5||$ll//-3/2(5/) + II $1155-3/2(90) + 


<C(u)[5 II $ II 75-3/2(5/) + h||$||^-3/2(5/) + S II $ II55—3/2(90) + 8‘ 


<C(u) (5 || $1175-3/2(5/) + ^11 $ II 55-3/2(5/) + 3 ' 
which in turn implies by taking S G M+ sufficiently small that 

ll$ll//-3/2(5/) < C(u)8. 


(4.51) 


Finally, by inserting (4.50) and (4.51) into (4.48), we immediately have (4.5). 


The proof is complete. 


□ 


Remark 4.4. For our study on the regularized full-cloak in the present section, 
we have made use of D$ = D$ U Dg U D\ as the virtual domain for the blowup 
construction. However, we would like to remark that one can also simply use 
Dg as the virtual domain by excluding the two “caps”, Dg and Dg. In such a 
case, one could replace Sg and Sg, respectively, by 5?&{Pq) and ^s(Qo) in our 
earlier arguments. Then by following a similar argument, one can arrive at the 
same conclusion as earlier for the approximate full-cloak. The reason that we have 
chosen to work with Dg instead of Dg as the virtual domain are two-folded. First, 
we would like to include the geometry considered in 28 for the regularized full- 


cloaks as a particular case in our this study; see also the discussion made after 


(4.8). Second, the corresponding mathematical arguments of dealing with D$ are 
more general than those of dealing with Dg, and we appeal to presenting a more 
general mathematical study. 


5. Regularized partial-cloak 

In this section, we consider the regularized partial-cloak by taking the generating 
set To to be an open subset on a flat plane Po in M 3 . Without loss of generality, 
we assume that Po is the {£3 = 0}-plane. If To is regarded as a domain in M 2 , 
it is assumed that To is bounded and simply connected with a convex boundary. 
However, in order to ease the exposition and illustration, we shall confine ourselves 
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to a special case by taking T 0 to be a square throughout the rest of the paper, 
which was actually considered in 28 . Nevertheless, at this point, we would like to 


emphasize that our subsequent study on the partial cloaking can be easily adapted 
to deal with a more general generating set r 0 . 

Let n 6 S 2 be the unit normal vector to r 0 and let q € M + . We next introduce 
the virtual domain D$ for the blowup construction of the regularized partial-cloak; 
see Fig. [2] for a schematic illustration. Let 



Figure 2. Schematic illustration of the domain D q for the regular¬ 
ized partial-cloak. 

D° r 0 (x) x [x — r • n, x + r ■ n], x € r 0 , 0 < r < q, (5.1) 

where we identify To with its parametric representation Fo(x). We denote by D v q 
the union of the four side half-cylinders and D q the union of four corner quarter- 
balls in Fig. [2j Moreover, we set S q and S q to be 

S] :=dD q n dD \, S 2 q : = dD q n dD 2 q . 

The upper and lower surfaces of D q in Fig. [2] are, respectively, denoted by 

rj := {x + q • n; x G r 0 }, T 2 q := {x - q • n; x G r 0 }. 

Define S° := Tj U Y 2 q . We then have 

D q = D° q U D\ U D\ and dD q = S° q U S\ U S 2 q . (5.2) 

Let 5 G M+ be the asymptotically small regularization parameter. We let D$ be 
the virtual domain and its boundary is clearly given by 

dDs = S° s U S] U ^ 2 . (5.3) 

In what follows, D$ shall be referred to as a screen-like region. If q = 1, we shall 
drop the dependence on q of D qi S q , 5)j and S' 2 , and simply write them as D,S°, 
S', and S 2 . 
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We shall make use of the blowup transformation A introduced in 28 that maps 
Ds to D , and refer to that work for the detailed construction. For our subsequent 
use, it is stressed that in D° s the blow-up transformation takes the following form 


A(y)=y ■= 


/ T 

f e 3 e 3 


I ^“ + e i e i + e 2 e 2 1 y> y e D s , 

where y E D$ and y E D, with the three Euclidean unit vectors given as follows 
e x = (1, 0, Of, e 2 = (0,1, Of, e 3 = (0, 0, if. 


Now we introduce the lossy layer for our partial-cloaking device as follows 

as(x) = <ti(x) := yB~ 2 (x), q s (x) = qi(x) := a + i/3 for xED s \D 5 / 2 , (5.4) 

where B(x) := V x A(x) is the Jacobian matrix of the blowup transformation A. 
The following theorem quantifies our partial-cloaking construction. 


Theorem 5.1. Let Ds be a screen-like region described above with its boundary 
giv en in (5.3). Let us be the solution to (2.13 ) corresponding to as and qs defined 
in (2.12), 


with ( Ds\D S / 2 ]cri,qi ) given by \5.Jj]) , and (D S / 2 -, cr a , q, 

'■* a ' ,x ' d be the incident plane wave satisfying 


being arbitrary 


but regular. Let u l {x) = e 1 ' 


|d ■ n| < e, e 1, (5.5) 

and let u^ix., d) be the scattering amplitude of us- Then there exists <5 0 E M + such 
that when 5 < <5 0 > ff^d) satisfies 

|i4,(x,d)| < C(w)(e + <J), (5.6) 


where C{oj ) is a positive constant depending on c o and D, but independent of a a , 
q a and x, d. 


Remark 5.1. Similar to Remark 4.1, one can immediately infer by Theorem 5.1 
that the push-forwarded structure 


{fl e ; cr, q} = a s , qs} 


with Fs dehned by (2.11) and {Ll e -,as,qs} defined in the theorem, 
approximate partial cloaking device. The incident aperture of the 
partial-cloak is given by 


produces an 
approximate 


Aj = {d E § 2 ; |d ■ n| < e}, 


whereas the observation aperture is A 0 = S 2 . 
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5.1. Asymptotic expansions. Similar to our notational usage in Section [4| if 
we let z denote the space variable on T 0 , then for any y G dD$, we define z y to be 
the projection of y onto r 0 . Denote by V n the normal derivative with respect to 
r 0 . Then for x G M 3 \ D$ we have 


S D s [ < />](^) = Gv;(x — y)fi(y)da y 
JdD s 

= / G u (x- Zy)4>(y)d<jy + 5 / G^- z y )fi(y) da y 


i s o 


's 1 


y 

21_l—l | 


(5.7) 


+ <5 J (y - z») ' V„,G„(x - 7.y)<p{y)duy + Op |x| \\<t>\\ H -yy eD )J, 

where and also in what follows, 0(y) := fi(y) and y := A(y) G D for y G Dg. 

We next estimate </>(y) in ( |5.7 ), and first derive the following lemma. 

Lemma 5.1. Let D$ be a screen-like region described earlier. Then for x G dD§ 
we have 

(1C" )*\(h \ixl = / (^5°)*[^K X ) + ^^W^Wh-o/^od))-, __ x G 5° U S 1 , 

[ D.) m i \ ( ^ o) .p ](i) + (Ksl) .p ](i) + 0(i |P|| H _ s/1(8D))! 


9“ UUxI = / 5 s»[0](*) + 0(i5|pilij-»/» ( 8£>))._ X e S° u S‘, 

1 ^so[^](x) + *Ss2p](x) + 0(511^11^-3/2(88,,), X £ S 2 , 


(5.8) 


where (/C^ o )*[0](x) is defined by 

(VM*M(X) -fif '0(y)* S , 

4:71 iso \ |Zj. Zy\ \Z X Zy | / 


and 


5^0 [0] (x) = / G u (z x ~ z y )fi(y)da y . 
J so 


Proof. By following a similar argument to that in the proof of Lemma 4.2 
compute for x G S$ U S} 


we 


|x y| = |x z x (y Zy) + z x - z y \ = \z x - z y \ + £>(<5), 


and 


(x - y, v x ) = (z x - z y , v x ) + 0(6), = e i ^~^ + 0(6). 

Note that for x G S'! one has 

|x-y| = 5|x-y|, e H x ~y| = i + 0(6). 

The proof can then be completed by using the similar expansion method as that 
in the proof of Lemma |4.2| □ 
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Next we introduce an auxiliary scattering problem for our partial-cloaking study. 
Let v$ be the solution to the following scattering system 

f A vs + u j 2 vs = 0 in M 3 \ D$ 
dv 5 


= 0 


on dDx 


dv 

u + := Vs — u 1 satishes the radiation condition (2.3) 


(5.9) 


Then we have 


Theorem 5.2. Let us be the solution to (5.9). Then there holds the following 
far-field expansion 


(vs -tt*)(x)= / G w (x - z s )K[n ■ Wu l (z)\(zy) da Zjj + 0(8), 
J r 0 


where 


K ; = ( \ ~ (*¥„)* (7+ (*¥„)' 


with the operator (K,f )* defined by 


(^fo)’M ; = P- v - / n ' V a G w (z — y)(j>(y) da z 
J r 0 


(5.10) 

(5.11) 

(5.12) 


Proof. The solution to (5.9) can be represented by (3.6), where the potential <f> 
satishes 

v 2 + W*)M = ~ on dD *- ( 5 ‘ 13 ) 


Using the fact that 


dv 


- + (££>)* : H~ 3/2 (S °) -> H~ 3/2 (S°) 


is 


invertible (cf. 1110 ). and by ( |5.13[ ) one has 

du l 


4 >( y ) — ~ (2 + (^s°) 


dv 


Next we can expand u 1 with respect to z as follows 

U*( y) = U l (z y ) + 8Vu l (z y ) • (y - z y) + C>(<5 2 ), y G dD s 


y G dD. 

(5.14) 

y e dD s . 

(5.15) 

y e r 1 , 

(5.16) 


By ( 5.14| ) and (5.15), one further has 

\ + (*%>)*) ['i ■ Vu‘(z*)] (y) + 0(6 

and hence ||0|| h- 3 / 2 (od) A C(ui). Noting that for y G T 1 , (y — 2n) G T 2 , we dehne 
V’(y) := 0(y — 2n) for y G T 1 and iffy) <f>(y + 2n) for y G T 2 . By using the fact 
that 

Vy- 2 ny = -Vy = -n, for y G T 1 , 
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one obtains (assume for a while that x G S °) 


4>(y -2n) = -0(y) = (j- - (/C^o)*) [n- Vu*(z £ )l (y) + 0(5), y G T 1 . (5.17) 


By inserting (5.16) and (5.17) into (5.7), there holds for x G M 3 \ D 
S£J<AK X ) = / GU X - Zy)(j){y)day + 0(6) 


i s o 


= / Gd(x - Zy)4>(y)d<jy + / Gd(x - Zy)i>(y)d(jy + 0(5) 


'r 1 


'r 1 


= J r gj * - * i ) ((j - (^»)-) _1 - (1 + ( K &) ' 

+ OW= [ G„(x-z s )K[a„ti , (z)](z s )<iCT a ,+0( l 5), 
«tr 0 

where d n u L (z) = n ■ Vu l (z) and K is defined by 

K := 2 (1 - (£$,)*) (Kg,)* (1 + (AC&) 


a n M*(z) (Zy)dcr z . 


-l 


(5.18) 


ft is noted in (5.10) that the function 0(z) := n • Vrd(z) depends only on z G r 0 , 
one can readily verify that 

(K? 0 )*[0](z) = l(A^,)*[9(y)](z + n), z e r„, 

for 0(y) = n ■ Vn*(z^), y G Id U T 2 . Hence, we can replace the operator K in 
(5.18) to be 

K := H - (K“ 0 )*) (AC?.)* (( + (K?„)* 


-1 


The proof is complete. 


□ 


Theorem 5.3. Let v§ be the solution to (5.9) and ^(x, d) be its scattering am¬ 
plitude. Then there holds 


4>( x , d ) + ^ 


-iouf E y i m (x)Y] m (z)|z| 


K[n • Vw l ](z) da z <C(uj)6, (5.19) 


where K is defined in (5.11). 


Proof. We recall the following addition formula (see, e.g., (37 


oo n 


l x — y| 


47r E E 


Jn 


2n + 1 




rn- 1-1 ' 


(5.20) 


n =0 m=—n 

where (r, £, d) and (r 1 , d') are the spherical coordinates of x and y, respectively; 

and Yff is the spherical harmonic function of degree n and order m. For simplicity, 
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the parameters (r,£,$) and (r,£ / ,'d / ) shall be replaced by (|x|,x) and (|y|,y), 
respectively. It then follows from ( |5.20[ ) that 

a 1 II 

l x -y| = l x l(i - y E n m (i)W)o) + ©(IxO). (5.21) 

m—— 1 1 1 

By (pl|), we have the following expansion 

(v s - u t )(x) = / G w (x - y)(/>(y)da y 
JdD s 

1 e^ |x| /■ E E-(x)y™(y)|y| 


(5.22) 


47T |x| 


'dD s 


0(y) + 0 (|x 


—2\ 


Expanding the integral kernel function in (5.22) around z y , one further has 

'</>(y )d(T y 




1 

47T 


-*•»f E I7"(*)yr(y)ly| 




E u 1 m (x)y 1 m (z s )|z 5 | 


'riur 2 


4>(y)day + 0(5), 


which together with (5.16) readily yields 

If -iuf E z 

^(x,d) = -— / e 


K[n ■ Vn*](z) da~ + (9(<5). 


'r 0 


The proof is complete. 


□ 


In view of Theorem 5.3, one immediately has the following result 


Proposition 5.1. Let vs be the solution to (5.9) and ^(x, d) be its scattering 
amplitude. Suppose that the impinging direction d of the corresponding incident 
plane wave u l satisfies 

|d ■ n| < e, e 1. (5.23) 

Then for sufficient small 5 € M + there holds 

|uoo(x, d)| < C(u)(e + <5), (5.24) 

where C(u ) depends only on T 0 and u. 


Proof. Let K be defined in (5.11). Since 

l|K[n ■ Vu‘(z)]|| H -i /2(ro) < C(w)||n ■ W(z)|| H -i /a(ro) , 
which together with (5.19) implies that 

boo(x,d)| < C(o;)(||n- Vu*(z)|| H -i /a(ro) + <5). 

The proof can be completed by noting that 

n ■ W(z) = ium ■ de iujz A . 
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□ 


5.2. Partial cloak. 


Proo f of Theorem 5.1. We follow similar arguments to those in the proof of The¬ 


orem 


4.1 


Define </>(y) := 0(y). Then it follows from 
I _ , r dus 


( 5 +ra,)*)i 


dv 


du l 

- — on dD s , 
dv 


and Lemma 5.1, and the proof of Proposition |5.1 that 

I|0||r-3/2(sj) < ||<l||^- 3 / 2 ( 5 i) + \\v ■ V«( z )|| H-z/ 2 (dD) + 3 — 0,1,2. (5.25) 

By® , there holds 


\u 6 -u l \ = \S Ds [(f)} \ = 


Is 0 


G«j(x - Zy)(j)(y)d(jy 


+£ ) ( ( KII011 n-3/2(so u s i )+^|| 011 H- 3/ 2 (ao))), 


which together with (5.25) further implies 


Woo(x,d)| < C(u)( || < h||^- 3 /2( S o)+e + 5 + 5||<h|| // ^3/2( 5 i)+ 5 H^Htf-s/a^)], (5.26) 


where $(x) = 4>(x) := ^ 

obtained by using a comple 
43 


(x) for x G dDs- The following estimates can be 

+ 

ely similar argument as that in the proof of Lemma 


m 


H-3/2(Si) <C(u)5 1 )/ 2 ||ll5||L2(D 5 \Di/2)’ J — 0,1, 2. (5.27) 


Inserting (5.27) back into (5.26), we have 




, d)| < C(uj) (5 1/2 \\us\\l2(Ds\d s/2 ) + e + ^]. 


(5.28) 


Following a similar argument to that in the proof of Lemma |4.4[ one can show 
(4.43). By change of variables in integrals, one further has 

f du t 5 


'dDs 


dv 


u l < 


Iso 


u*(z)$(y )d(Jy + C(u)5\\$\\ H -3/2 {dD) 


(5.29) 


<CM||$|| ff -3/2( 5 0) + C (Uj)5\\$\\ H -3/2(Q D ^ 


which together with (5.8) readily yields 


I dD s 


dv} 

-^{u 5 -vf)\ + 


I dDs 


r)iii _ 


< 


du 


——( z )55o[0] + (7(5(11$ || H-3/ 2 {dD) + 1) 


(5.30) 


Iso dv 


< C(u)e + C'^(|| < h||//- 3 / 2 (aD) + !)• 
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By taking the imaginary parts of both sides of (4.43), and using (5.27), (5.28), 
(5.29) and (5.30), there holds 

W U s\\l^(D s \D s/2 ) 

<C(o;)|u 00 | 2 + C(lu) 11 11 #-3/2(50) + C(u>)e + C(o;)5(||l ) || H -3/2( 9D ) + 1) (5.31) 

<C{lo) ( '8\\u6\\l2( D5 \d 6/2 ) + ^ 1/2 \\ u s\\l 2 (Ds\d s/2 ) + e + 

Then by choosing a sufficiently small 8 e R+, one easily has from (5.31) that 

\\us\\l^(d s \d s/2 ) < C(co)(e 1/2 + 8 1/2 ). (5.32) 

Plugging (5.32) into (5.27), one finally has 

II^ > IIr- 3 / 2 (5°) < C{u)(e + an d II^IU-VRao) < C{ui){e ^ 2 + <X^ 2 )> 

which together with (5.26) readily yields (5.6). 

The proof is complete. □ 
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